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Abstract. We consider the general Z2-symnietric free-fermion model on the finite 
periodic lattice, which includes as special cases the Ising model on the square and 
triangular lattices and Z„-symmetric BBS r'^'-model with n ~ 2. Translating 
Kaufman's fermionic approach to diagonalization of Ising-like transfer matrices into 
the language of Grassmann integrals, we determine the transfer matrix eigenvectors 
and observe that they coincide with the eigenvectors of a square lattice Ising transfer 
matrix. This allows to find exact finite-lattice form factors of spin operators for the 
statistical model and the associated finite-length quantum chains, of which the most 
general is equivalent to the XY chain in a transverse field. 



PACS numbers: 75.10Jm, 75.10.Pq, 05.50+q, 02.30Ik 
1. Introduction 

Consider a two-dimensional M x N square lattice with spins cr = ±1 living at each of 
its sites. The most general Z2-symmetric plaquette Boltzmann weight is given by 

W{ai, 0-2, as, (74) = ao(^l+ ^ aijCiaj + aAaia2(y'ia^ . (1) 

l<j<j<4 

Partition function of the corresponding statistical model can be represented as an 
integral over four-component lattice Grassmann field with a quartic interaction [6l [7] , 
which disappears if the parameters satisfy the condition 

04 = 012034-013024 + 014023- (2) 

The model ([1]) can be mapped to the eight- vertex model in an external field [2]. The 
analog of the free-fermion condition (E]) in the vertex picture was obtained earlier in 

Below we study the model defined by ([I])-(l2]). It appeared in several equivalent 
formulations in papers by different authors. Its partition function can be evaluated 
by a variety of methods for both infinite and finite lattice [U [6l [TOl [131 ISO]- Elliptic 
parametrization of the Boltzmann weights and inversion relations for the partition 
function were established in [3l [5]. Grassmann integral representations for the 
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corresponding transfer matrix and its eigenvectors have been found in [21]. Infinite- 
lattice correlations in a dual model were studied in [J^. The model (II])-(I2]) includes 
as special cases the Ising model on the square and triangular lattices, as well as Z„- 
symmetric Baxter-Bazhanov-Stroganov r'-^-'-model with n = 2. 

The transfer matrix 14[cr, c'] is a function of 2N spin variables ao, ■ ■ ■ , (Tn-i, 
o"q, . . . , cr^_i given by the product of plaquette weights ([1]) over one lattice row (Fig. 1): 

N~l 

VeW,a'] = H W{aj,a'^,cr'^^„a,+,). (3) 

j=0 

The subscript e = ±1 corresponds to periodic and antiperiodic boundary conditions on 
spin variables in the horizontal direction, i.e. = eao and cr^ = eo-'q. 

-.^1 -.-I 



G, Qi G2 O-j Qj+i G^_2 a^.i 
Fig. 1 

The transfer matrix acts on the 2^-dimensional space of maps / : (^2)^^ — t- C by 
{VJ)[a] = J2Vs[cr,a']f[a']. 

It commutes with the operators of translation and spin reflection defined by 

iTef){cro, . . .,aN-i) = fieri, ■ ■ . , (Tat-i, ecro), 
{Uf)[a]=f[-a]. 

Partition function of the model can be written as Z = Tr ^V^^'^f/^~2^ j , where 
e' = ±1 corresponds to periodic and antiperiodic boundary conditions in the vertical 
direction. To write spin correlation functions, one should introduce spin operators {sj} 
with J = 0, . . . , — 1 defined by (sjf) [a] = o-jf[a]. Then for ji < ^2 < • • • < jn we have 

where 

s^,k = ViskV-^ = VJT^soT-'V~^. 

The computation of correlation functions therefore reduces to finding matrix elements 
(form factors) of the spin operator sq in the common basis of eigenstates of V^, and 
U. The present work is devoted to the solution of this problem. 

Multiplication of the Boltzmann weight ([1]) by e^^'^i°^2-cr3cr4) ^jth any K does not 
change the transfer matrix. Thus, up to overall factor, nontrivially depends on five 
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parameters. The parameter set may be thought of as a five-dimensional projective space 



with homogeneous coordinates 

K = (ai2 + 034) («i3 + 024) + (ai4 + 023) («4 + 1), (4) 

A = (ai4 - 023) (04 - 1) - (ai2 - 034) (ai3 - 024), (5) 

= (04 + if - (ai2 + 034)^ - (ai3 - 024)^ + (ai4 - (^23)"^, (6) 

g = 4 (014023 - 013024) , (7) 

r = (04 + 1)^ + (012 + 034)^ + (oi3 + 024)^ + (ai4 + 023)^ (8) 

V = (O12 + 034)(Oi3 + O24) - (ai4 + «23)(«4 + !)• (9) 



It will be shown below that this set is foliated by two-dimensional surfaces such that the 
transfer matrices corresponding to different surface points are mutually commuting and 
diagonalized in the same basis. Different surfaces are labeled by the triples ^, 
Each of them contains a curve representing a set of transfer matrices of the Ising model 
on the triangular lattice, and this curve contains a point corresponding to the transfer 
matrix of an anisotropic Ising model on the square lattice. 

As a consequence, finite-lattice form factors of the free-fermion model ([I])-© can be 
obtained by a suitable parametrization from the Ising ones, for which the corresponding 
formulas were conjectured in and proved in [II1II2]- A more elegant proof based on 
earlier results of [HI |23] and Frobenius determinant formula for elliptic Cauchy matrices 
was recently found in [T6] . 

It is well-known that the six- and eight-vertex model transfer matrices commute 
with the hamiltonians of the quantum XXZ and XYZ spin-1/2 chains, respectively. 
Similarly, the diagonal-to-diagonal transfer matrix of the Ising model on a square lattice 
commutes with the hamiltonian of the quantum Ising chain [22]. Local hamiltonians 
are also known for the Ising model on the triangular and hexagonal lattices [25] . 

We will show that to every above-described two-parameter family of commuting 
transfer matrices of the general free-fermion model one can associate a local spin chain 
hamiltonian depending on three parameters ^, ^, ^. It turns out to be related to the 
hamiltonian of the quantum XY chain in a transverse field by a similarity tranformation. 
Simple considerations then allow us to rederive spin form factors of the finite-length XY 
chain, recently obtained in [15] by the method of separation of variables. 



2. Kaufman's approach and Grassmann integrals 

Let us define two operators Vaphj the relations 

Since V± and Va^p commute with U, the set of eigenstates of V+ (V-) consists of even 
(odd) under spin refiection eigenvectors of Va and odd (resp. even) eigenvectors of Vp. 

It was observed in [21] that VJj.p can be naturally represented as 2A^-fold Grassmann 
integrals. To describe this result in more detail, introduce auxiliary Grassmann variables 
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ipo, ... ,11)^-1, ipo, . . .ipN^i. Our convention will be that "dotted" variables commute 
with the usual ones and anticommute between themselves. We will also need to consider 
discrete Fourier transforms of variables of both types 



J=0 



using two sets of quasimomenta: Oa = {^,^,...,27r — and dp 
{O, . . . , 27r — In this notation, one has [21] 



N-l N-1 



V^[a, a'] = Cu [ VtpVijj JJ e"^^^ JJ e"^'''% = a,p, (10) 



j=0 j=0 

where = [Ue&e, Xe] ' and 



l2 



Xe = [ai2 + ^34 + (ai3 + «24) cos 6f + [(aia - 024)^ + 4014023] sin^ 6, (12) 

X6»G'ii(^) = 22 sin 6 [023 + 012024 + 013034 + 014O4 - 2 (014023 - 013O24) cos 6] , (13) 

XeG22{0) = 2i sin 6 [014 + O12O13 + O24O34 + O23O4 - 2 (014O23 - O13O24) cos 9] , (14) 
XeGi2{9) = XeG2i{-9) = [(012 + 034) (04 + 1) - (013 + 024) (ou + 023)] 

+ [(0i3 + 024)(04 + 1) - (O12 + 034)(0i4 + O23)] COS 6* (15) 

+ [(024 - 013) (04 - 1) + (012 - 034) (oi4 - 023)] z sin 6*. 

The representation f lTU]) was used in |5T] to find the eigenvectors of Va^p explicitly 
in terms of spin variables. One of the reasons to proceed in this way was the difficulties 
with the well-known Kaufman's algebraic approach [12] designed for the two-dimensional 
Ising model. It can be expected on general grounds that the conjugation by Va^p induces 
linear transformations of the standard Clifford algebra generators. However, in contrast 
to the Ising case, it was not clear how one can determine the explicit form of these 
transformations. 

This can be done as follows. First look for a representation for the inverse V~^[a, a'] 
in the form f lTO|) : 

/7V-1 N~l 
e ^-[^'^1 Yl e"^'^^- Y[ e"^'^^' , (16) 
j=Q j=0 

where the overall coefficient d, and quadratic action Siy[ip,ip] should be chosen so that 

W'] j=o 
Using ffTOl) and the ansatz f|T6|) . the l.h.s. of the last relation can be written as 

/N-l N-l N-l N-l 

la'j i=0 j=0 j=0 j=0 
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Let us explain how this expression can be simpHfied, since below similar calculations 
will be done without going into the details: 

• Only even part of 11 jl^^ e'^^^^ H^LV e"^^^ gives non-zero contribution to the integral. 
Since it commutes with any quadratic expression in Grassmann variables of both 
types, the factor e'^"''^''^! can be put together with e'^"''^''^!. 

• The summation over intermediate spins can then be easily done: 

N-l N-1 



[a'] i=0 i=0 

where we use shorthand notation ij)(p = J2j^=o '^j'^j- '^'^^^ transforms (fT8|) into 

/N-l N-l 



j=0 j=Q 

To pull e^'^ to the left through = nj=o^ e'^^^^ , write the latter factor as a sum 
of its odd and even part, which gives F[^]e'^'^ = e'^'^Fevcn[^] + e~'^"^Fodd['^]- In the 
integral corresponding to the first term, make the change of variables "0 — )■ —i/j, 
ijj — )■ —ijj. Since Slip, ip] is even, this allows to rewrite the previous expression as 

7V-1 N-l 



^ A — n A — n 



3=0 j=0 

Thus, in order to satisfy ( fTTl) . it is sufficient to choose and Sy[ip, ip] so that 



In view of (fTT]) . it is natural to try to satisfy this relation with Siy[ip, ip] diagonalized by 
discrete Fourier transform. This indeed works and one finds that 



det G{e) 

where the 2 x 2 matrix G{e) is defined by f in]) -f lT^ . 

The next step is to introduce the Clifford algebra generators {jOj}, {qj} {j 
0,...,iV-l), 

\ i-i N-l / I 

which satisfy standard anticommutation relations {pj-,Pk} = {QjiQu} = '^^jk, {Pj^Qk} = 
0. It is easy to verify that these operators can be represented in a form similar to 
Grassmann integral representation ( |TOl) : 

> [a, a'] = 2-^ / Vr]Vfj e"'^ ."^^ J] e"^^"^' J] e^^K (19) 

J \ i[rjj r]j) I 



k=0 k=0 
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We especially note that Fourier transforms pg, qg are obtained by simple replacement of 
the subscript j by 6 in the r.h.s. of f|T9|) . 

It then becomes possible to compute the result of conjugation of {pj}, {qj} by V,y 
using the representations (fTOj) . (fT6l) and (1191) . After summation over intermediate spin 
variables one obtains 



iiVe - Ve) 



j=0 j=0 



The integration over ip, rj, r], ip can be performed in the Fourier basis. Comparing the 
resulting integral over i/j and ip with (fT9l) . we find that the induced rotation is explicitly 
given by 

(20) 

ae = X9[l- det G(^) + G'n(^) - Gssl^)] , 
/3, = [1 + det Gl^) + Gn(^) + GsaC'^^ 



where 



One can now follow Kaufman's method and find the eigenstates of p by diagonalization 
of the two-dimensional rotations fl?Il) . 

It is straightforward to check using (fT2|) - (fT5|) that ag, j3g can be written in terms 
of the parameters (jl])-© as follows: 

ae = T + 2vcose + 2iX sin 9, /3g = -ge^'^ + 2Ke'^ - /x. (23) 

and that aga-e — PgjS^g = 4x^6*12 (6') 6*21 (6*) • In fact the relations fj23|l are the origin of 
the parametrization (j4])-(l9l). Let us further introduce three parameters /Co, JCx, by 
A sinh 2/Co At cosh 2/Cj, + cosh 2/Co Q cosh 21Cy — cosh 2/Co 

K cosh 21Cx sinh 2/Cj, ' k cosh 2/Ca; sinh 2/Cj^ ' k cosh 2/C2. sinh 2/Cy 
It will be shown in the next section that the eigenvectors of Ve coincide with the 
eigenvectors of a non-symmetrized transfer matrix of the anisotropic Ising model on 
the square lattice: 



exp 



E exp{K:. E C,] exp{^!^ (24) 



2 

i=0 i=0 i=0 



where tanhA^* = e and the operators {Cj} are defined by 

{Cjf){- ■ . ,aj^i,aj,(Tj+i, ...) = /(.. . ,(Tj_i, -crj,aj+i, . . .). 

We assume for concreteness that JCq, IC^, fCy are real and positive and /C* < fCy. This 
mimics the ferromagnetic region of Ising parameters. 
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Define the functions 

satisfying bfb^g = 1. The square roots here and below are fixed by the requirement of 
positivity of the real part. Determining the eigenvectors of the induced rotations ( |2T]) . 
introduce the creation-annihilation operators 




-ie 



'e P-e 



-i\ b^ 



•Q- 



Pe 



AO 



+ iy 



Qe 



(25) 



where 



P-0 



V2 



btb" 



bab' 



These operators satisfy canonical anticommutation relations {'0^,^^/} 
{ipe^'^B'} = 0, {^^,-061/} = defii. They transform diagonally under conjugation by Vy. 





eee„ 



where 



8o = In 



4xeGi2(^) 

The matrices Vv can therefore be written as 



K = 2^a^[n^^^i2(^) 



exp 



(26) 



(27) 



where the overall scalar multiple can be fixed e.g. by the identification of eigenvalues 
with formulas in [21]. In contrast to the Ising case, Eq ^ £^q and hence the spectrum 
is (generically) nondegenerate. In particular, the eigenstates of will automatically 
diagonalize the translation operator. 

The left and right eigenvectors of are multiparticle Fock states 



{21 



^ {vac\ipg^ . ..ipB^, \9i, . . .,9k)u = ■ ..tpljvac)^, 

where 0i, . . . ,6^ G 6i,- The corresponding eigenvalue is equal to exp | ^ Xleee^ ^e ~ Si=i ^f; 

Vacuum vectors \vac)iy and u{vac\ are annihilated by all i/jg (resp. ipl) with 9 E Oi, and 
are normalized as ^{vac\vac)i, = 1. Since V is not symmetric, i,{vac\ is not necessarily 
hermitian conjugate of \vac),y. 

Not all of the states (l28l) are eigenvectors of the full transfer matrix (and the 
associated quantum spin chain hamiltonians below). For periodic and antiperiodic 
boundary conditions on spin variables the number of particles in these states should 
be even and odd, respectively, and one has 



^1) • • • ) ^2fc)a, 



T+I^l, . . . , 92k) a,p — 

T-\9i, . . . , 92k+l)a,p = e~ ^^=1 '\9i, . . . , 92k+l)a,p- 



2fc + l , 



(29) 
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We especially note that the transfer matrix eigenvectors and form factors depend 
only on /C^., Ky, ICq and are independent of ^, ^. The latter two variables appear only 
in the eigenvalues and can be thought of as spectral parameters. 

3. Form factors 

Consider instead of Vs the conjugated transfer matrix = SV^S"^ with S = 
exp 1^ 'l2f=o "^^^ matrix of induced rotation in fl20|) then becomes A{Vl) = 

A-\S)A{V^)A{S) with 

cosh /Co isinh/Coe~*^ 
-i sinh /Co e*^ cosh /Co 



AiS) 



and it can be straightforwardly checked that A(yi) = A{Viy) . Therefore the 

^0=0 

orthonormal system of left and right eigenvectors of can be chosen as follows: 

■ ■ ■ ,dk)u = S~^\\9l, . . . ,6k)u) , ui^l^ ■ ■ ■ i(^k\ = ■ ■ ■ ,dk\) S. 

V / Ko=0 V //Co=0 

These vectors are of course proportional to those in ( l28l) . 

Since S commutes with the operators {sj}, nontrivial spin form factors coincide 
with those computed for /Cq = 0: 

J^m]n{^^^') = f(^l) • • • ,(^m\si\6[, ■ ■ ■ , 9n)p = [p{0'i, . . . ,9'Jsi\9i, . . . ,9m)a] = 

= [a{9i,...,9m\si\6[,...,6'^)p]j^^^^Q, (30) 

where m, n are simultaneously even or odd. On the other hand, setting /Co = leads to 
important simplifications in the formulas of the previous section. It implies in particular 
that X = 0, ae = a^g, = bg and pe = 1. Moreover, 



1 - tanh /C* coth ICyC'^) (1 - tanh /C* tanh ICyC''^) 
1 - tanh /C* tanh /C^e*^) (1 - tanh /C* coth /Cj^e"^^) ' 

and thus the Fock states (l28l) with /Co = coincide with the eigenvectors of the 
symmetric transfer matrix of the anisotropic Ising model on the square lattice, cf. 
[16]. Let us therefore denote \9i, . . . , 6'fc)[f™^ = . . . , 6*^)1/^ • In the symmetric 
Ising case, the left and right eigenvectors are related by hermitian conjugation, i.e. 

= {\9,,...,9kyr^y. 

To write the corresponding form factors, introduce two functions 7e, ue and two 
parameters ^, defined by 

cosh 75) = cosh 2/C* cosh 2/Cy — sinh 2K.I. sinh 2/Cy cos 9, > 0, 

n.,^. sinh^^^ 
ue = In . , e = [1 - (sinh 2/C,. sinh 2]Cy)-'] * , 

Lie'eep ^^^^ 



2 



n n ' 



2 7e+7fl' 
2 



ne,e'e..sinhl^ne,e'g..sinh2^ 
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-i(l-\)e,+uoj2 _n_ m~\)B',^Ve>n 



^gj^,^o=^^— "^v^n ■ / n / X (32) 



X 



/ sinh 2K,y \ " tt sm ^ sm -p-^ smh — 



. , 70^+79,- 11 7e'+7£)'. 11 



\sinh2A^ / • J-A , 7e'+7e'. 11 . ' 

^ l<j<j<m Sllin 2 l<i<i<n smh — — - \<i<m,\<i<n Sm — ^ — 

As explained in the Introduction, this formula allows to compute (long-distance 
expansions of) an arbitrary spin correlation function in the statistical model ([I])-(j2]). 

4. Quantum spin chain hamiltonian 

The matrix Vy clearly commutes (cf (1271) ) with the operator 

which has the following form in terms of {pe}-, {<ie}'- 

Hy = -^'^^{a0- a_e) {P-ePe - q-eqe) - ^iPe Q-ePe^- 

Now using fl22|) and performing the inverse Fourier transform, one finds that 

N-\ 

Hu = - (pjPj+i - qjqj+i) + ifiqjPj - 2iKqjPj+i + z^g^pj+a}, 

j=0 

where Pj+N = Pj, qj+N = qj ior u = p and Pj+N = —Pj, qj+N = —qj for u = a. 
Using the standard realization of {pj}, {qj} in terms of Pauli matrices cr^'^'^ 

Pj = (g) ■ ■ ■ (gXT^'^ Oer^ (g) 1 (g) . . . (g) 1 , 

j times N—l—j times 

qj = (g) . . . (g) a"" (g (g 1 (g) . . . (g 1. 



V 

j times N~l~j times 



one finds the quantum spin chain hamiltonian, which commutes with the transfer matrix 
Vs and the operator ( l24l) : 



N-1 



where the boundary conditions are o"J^_jv = crj, crj_,_^ = eaj, Uj+N = £<Jj. To put this 

1 i_ 

hamiltonian in a more familiar form, conjugate it by the matrix S' = VxVyS with 

K = exp j/c: = exp {/c, ^1^1+1}, 

i=o j=o 
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which gives 

N~l 

H'^ = S'HeS'-^ = const ■ {e"^^^a|(T|+i + e^^^ajaj+i + 2 coth 2ICy aj} . (33) 

j=0 

One recognizes here the hamihonian of the XY chain in a transverse field. If we 

1 1. 

restrict ourselves to the symmetric Ising transfer matrix by setting = Vy^V^Vy^, 

this reproduces the well-known observation of [26] that H'^ commutes with the second 

1 1 

symmetric Ising transfer matrix Vx? VyVx ■ 

The orthonormal eigenvectors of H'^ can therefore be chosen as follows: 

1^1, ... , ^fc)^ = e'^'^vf^vf^ 1^1, ... , (34) 
...M= '"^'(^1, • • . ,4|n^V.^^e±^, (35) 



where 



e-^W = exp{i5^7.-E>J (36) 



1 1 



is the corresponding eigenvalue of YxVyVi . Note that different signs in (l34l) and 
fl35|) give equivalent representations of the same vector, and we have the hermitian 
conjugation identities {6k, • • • , ^i| = (|^i, . . . , 

The operators erf = ipiQi commute with the Z2-charge U and thus have non-zero 
form factors only between the states of the same type (a- or p-). Since they are bilinear in 
fermions, their matrix elements are easily computable. The only nontrivial form factors 
of the XY chain fl33l) are those of the operators erf and af , which map the a-sector to 
j9-sector and vice versa. Define 

Then, using different representations for the eigenvectors flM|) -f l35|) . one finds that 

±zsinh/C:j-^Y„(0,0Vf) + cosh/C:j-„^j;(0,0Vf) = e^'-"^ X 

X (^1, . . . ,^„,|Kf V.^^(±^sinh/C:af + coshlClat)vfwf^\9[, . . . , C)?'^^ = 

I 7(e)-7(e') T . +i zri , , T ■ I l{B)-l{0') n\ 

= e^^^ • • • , e^lVy^^afVy^^ \9[, . . . , eX'""' = e^^^J'^liO, 6'). 

From these relations one obtains the final finite-length XY form factor formulas: 

sinh^^tom 

^Vf) = . ■ ' ^^{0, 0% (37) 
z smh /C* 

cosh^^tom 

•^S(^'^Vn = cosh^'- J'lnUO^O'), (38) 

where J^^]n{e,e') is defined by §^ and 7(6>) by ([36D. The same expressions have been 
recently found in [15] by the method of separation of variables and used to rederive the 
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asymptotics of two-point correlation function in the disordered phase without the use 
of Toephtz determinants and Wiener- Hopf factorization method. 

Detailed discussion of the XY chain is beyond the scope of the present paper. 
We would like to mention, however, the works [T71 [18], where time- and temperature 
dependent two-point correlation functions of the finite XY chain were expressed in terms 
of 2N X 2N determinants. In fact equivalent representations can be simply deduced from 
the observation that spin operators, as well as the densities e"^^", are elements of the 
Clifford group; the entries of the corresponding determinants are then given by thermally 
dressed two-particle form factors. We will report on these issues in a future publication. 

Our last remark concerns the operators appearing in the hamiltonians and H'^. 
It can be checked that they obey Onsager algebra relations for the standard generators 

= 2kAi + ^iAq - 2XGi + pA2, 
H'^ = const ■ (e^'^^Ai + e~^'^^A^i + 2 coth 2/Cj,Ao) . 

The similarity transformations with exp(A^ylo), exp{]C'Ai) relating these two 
hamiltonians are combinations of elementary Onsager algebra automorphisms described 
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